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Abstract. Fermi gases with short-range interactions are ubiquitous in 
ultracold atomic systems. In the absence of spin-flipping processes the 
number of atoms in each spin species is conserved separately, and we 
discuss the associated Ward identities. For contact interactions the spin 
conductivity spectral function a s (oj) has universal power-law tails at 
high frequency. We derive the spin /-sum rule and show that it is not 
j^JTJ affected by these tails in d < 4 dimensions. Likewise the shear viscosity 

spectral function rj(uj) has universal tails; in contrast they modify the 
viscosity sum rule in a characteristic way. 

cd; 

0" 1 1 Introduction 

Condensed matter systems near a phase transition generally have universal low-energy 
properties, while the high-energy response depends on non-universal details of the 
microscopic interaction. Ultracold atoms provide an important exception: in dilute 
gases, where the range of the interaction r e is much shorter than the mean particle 
spacing, also the high-energy properties are universal up to a cutoff energy h 2 lmr\ set 
by the interaction range, which can be much larger than the Fermi or thermal energies 
PQ . The correlation functions have characteristic high-frequency and momentum tails 
which are controlled by the Tan contact density C [SIEISIESIE] • This quantity measures 
the probability of finding two atoms of different species near each other. Together, two 
atoms can absorb a large kinetic energy and undergo a boost in opposite directions 
while conserving total momentum. Hence, the high-energy response of the system is 
proportional to the density C of such pairs. 

In this work we look in particular at the response to a magnetic field gradient, 
the spin conductivity a s , which has recently been measured [7] and provides an ex- 
ample of quantum limited transport. Aspects of this are understood within kinetic 
theory [§1[TU]. while a recent strong-coupling Luttinger-Ward calculation [5] explains 
the spin diffusion quantitatively and predicts the full frequency dependence of the 
spin conductivity a s (uj). Furthermore, we consider the response to shear flow, the 
shear viscosity r](ui) [TTlll2[H~51H^] in two and three spatial dimensions. The transport 
coefficients exhibit universal power-law tails at high frequencies, and we study how 
these tails affect the exact sum rules which link the frequency integrated response to 
the thermodynamic properties of the system [51 [LH[To1HBlll7U18l[T^] . The question of 
spin transport is connected with the conservation of the particle numbers N a of the 
spin species. For the case of a density-density interaction there are no spin-flipping 
processes and each N„ is conserved separately. This implies spin-selective Ward iden- 
tities for every spin species, which we will then use to derive the spin sum rule. 
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2 Model and symmetries 

Consider a two-component Fermi gas with contact interaction which is described by 
the grand canonical Hamiltonian 

H= [ ^ ^ (5 T - ^ - Pa) + ffo^Vjl^t (!) 

with mass m, chemical potential for spin species a =t-i? an d 7i = 1. At low 
energies s-wave scattering is allowed only between opposite spins by the Pauli prin- 
ciple. The contact interaction go leads to ultraviolet divergences which need to be 
regularized PQ. 

The interacting Fermi gas <m is invariant under a U(l) x U(l) symmetry corre- 
sponding to the separate conservation of f and J, particle number. This is readily seen 
by coupling a different gauge field to each spin component [3D]. In the absence of a 
magnetic field the symmetry is enlarged to SU(2). The spin-selective particle number 
and current operators can be written as 

p a (x) = ^ a {x)^ a (x) (2) 
3i(x) = -fyUx)M*(*) ~ flfVtolMaO) (3) 
where x = (x, t) . These operators satisfy the continuity equation 

d tP „(x) + d l fAx) = d^(x) = o (4) 

separately for each spin component with conserved particle numbers Nf and iVj. . The 
bare current operator ([3]) acquires no interaction correction in the case of the density- 
density interaction (|TJ) since \p a (x), p a * (y)} = and [H- mt — p a N a , p a i (x)} = 
|21) . The continuity equation implies spin-selective Ward identities which connect the 
number p a and current f a response functions with the Green's functions. These have 
been derived by Behn [22] . 

d,(T^(x)^(y)^l(y')) = S aal (TMy)^(y')) [S(x-y)-6(x-y')} (5) 
or in momentum space 

ujp a (p,cr',e;p + q,a',€ + uj) - qij l a (p, a', e; p + q, a',e + uj) 

= <W [G a (p + q, e + u) - G„ (p, e)] (6) 

with Green's functions G~ 1 (p,cu) = —uj + e p — p a — E a {p,uj) and the free single- 
particle dispersion e p = p 2 /2m. In particular, there is no response of the I Green's 
function to a q^j^ perturbation, which is not immediately obvious from looking at the 
perturbative contributions: indeed, Maki- Thompson and Aslamazov-Larkin vertex 
corrections [12] cancel exactly in this case. For SU(2) invariant models there is an 
additional Ward identity for the a + operator [32] . 

3 Spin /-sum rule 

The correlation functions of number-current xj n and spin-current Xjs are defined as 

/>oo 

Xjn/js (q,a/) = i/ dtdxe^-^^Cx,*), £,.(0,0)]) (7) 
Jo 
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in terms of the number and spin current operators j l n / s {x) = j\{x) ± j\{x). The 
corresponding number and spin conductivities in the zero-momentum limit are defined 
in terms of the correlation function ([7]) as 

lmxjn/ js (0,^) 

A Kramers-Kronig transformation relates the frequency integral of a n / s (uj) to the 
current correlation function at zero frequency, 

/OG 
C n /sM = Xjn/js(0,W = 0), (9) 
-oo " 

which is real. The Kubo formula ([7]) can be expressed in terms of the fermionic Green's 
and response functions in the Matsubara formalism as [50,13 

X js (0,0) = --L £ ^ x r* CT ,j*(pX,ie n ) = -± £ |j*(p,^e n ) (10) 

where r^ CT ,p^/m is the bare spin-current response vertex with Pauli matrix t z . This is 
multiplied with j^(p, er', ie n ), the fully dressed current response function from Eq. (O 
in the limit w = 0,q-)0. The Ward identity ([6]) for each spin component relates the 
current response function in the static limit lo = 0, q — ¥ to the Green's function, 

■ z , . x dG a (p,ie n ) p z dG a (p,ie n ) 

3,<P,">«n) = = ~- dp ■ (11) 

The Matsubara sum over the Green's function yields the momentum distribution 
— J2e n ^(Pj * e «) = n p<? an d one obtains 

= -VT, ^ = «.(0,0). (12) 

per 

The same result is obtained if one considers not the spin-current but the number- 
current with bare reponse vertex 8 aa ip z lm. The normalized integral over the solid 
angle Qd yields J dflp 2 z j fid — P 2 /d in d > 2 dimensions. Integration by parts over p 
then gives 



m f Qdd PP d 1 p dn p° _ \ " n p° l f?d ^ d A 

Xj„/jsiu,uj- (2 ^ d md Qp -2^ m md{27T) d2^P n p- 





(13) 



where we have explicitly written the ultraviolet momentum cutoff A ~ 1 / 1 r- e | . The 
first term gives the density, while the second term depends on the cutoff. For zero- 
range interactions the momentum distribution at large momenta is proportional to 
the Tan contact density, n p(J — C/p 4 as p — > oo in any dimension 013113] • Hence, 
the cutoff term CA d ~~ 4 vanishes for A — > oo (r e — > 0) in any dimension d < 4. In 
combination with Eq. ((5J) this completes the derivation of the particle number and 
spin sum rule 

1 [°° n 

- / duja n/s (uj) = — (14) 
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with the total density n = n-f + nj,. 

In the Galileian invariant model (JlJ the number current is proportional to mo- 
mentum and cannot decay. This implies the conservation of the total number current, 

[H,jdxji(x)]=0, (15) 

and consequently the number conductivity has a sharp Drude peak at zero frequency, 

77 

a n (ui) = — ttS(uj). (16) 

771 

In contrast, the global spin current is not conserved because scattering transfers 
momentum between f and J, particles, 

[H,Jdxji(x)]?0, (17) 

and the spin conductivity a s (uj) has a finite and nontrivial response at uj > 0. The spin 
conductivity in 3d has recently been computed in the Luttinger-Ward formalism [5]: 
there is a broad Drude peak at low frequencies, followed by a universal high-frequency 
tail 

<t s (o;-kx)) = o 7 ~ \ 3 / 2 ( 3d ) ( 18 ) 

in accordance with results from the operator product expansion [T7]. Both in two 
and three dimensions the tail decays sufficiently fast for the frequency integral (I14|) 
to converge, so again the universal high-energy properties of the zero-range model do 
not affect the form of the spin /-sum rule in d < 4 dimensions. 



4 Shear viscosity sum rule 

The shear viscosity 77 measures the friction of a fluid subject to a shear flow of both 
spin species simultaneously (mass flow). The frequency dependent shear viscosity 

1( U ) = (19) 

UJ 

is defined via the stress correlation function 

Xxy*vi<hu)=i rftdxe^ t - q - x )<[7J iC1/ (x ) t),7J xy (0,0)]) (20) 
Jo 

at zero external momentum, q = 0. The stress tensor operator has the off-diagonal 
components [2TUT3] 

^T~ c p-q/2,<T C P+ciA ff ( 21 ) 

p 

since the interaction correction vanishes in the zero-range limit |21U13j . Again a 
Kramers-Kronig transformation relates the frequency integral of the viscosity to the 
stress correlation function at zero external frequency (static limit), 

/oo 
V(u) = Xxyxy(u = 0). (22) 
-00 7T 
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In analogy with the spin case, the Kubo formula (1201) is expressed in terms of the 
stress response function T xy as }13) 



Xxyxy{0) — py ^ ] m T xy (p : ie n ) . 



(23) 



In the static limit of external ui = 0,q — > the stress response is determined by the 
Ward identity associated with momentum current conservation 23 , 



~Px- 



d G(p, ie r , 

d P y 



p x p y dG(p 1 ie n ) 
p dp 



and hence 



pIpI dn 



w,(o)=4e 

V rap op 



pa 



(24) 



(25) 



The normalized integral over the solid angle yields J dfldp\p 2 y l = P 4 /[d(d + 2)] 
in d > 2 dimensions. Performing an integration by parts as in Eq. (|13p relates the 
correlation function to the kinetic energy density Skin = y H pCT e p n pu- The integrals 
are ultraviolet divergent and can be regularized by a momentum cutoff A, 



Xxyxy (0) 



-Ekin 



md(d 



•2) (2*) 
„4 



Through the momentum distribution n pa = C/p for p 
term depends on the contact density C, 



Xxyxy (Q) — t -Ski] 



CA 



d-2 



(2ir) d md(d + 2) 



(26) 

oo (see above) the cutoff 
(27) 



The kinetic energy density Sdn can be re-written using the Tan relations for the 
internal energy density e or the pressure P as [2,3 



Skin — £ 
Skin = £ 



c 

Aura 

C 
Aura 



ln^ = P- 

£B 



C 
Anm 



C UJ A 
In — 

Airm Eb 





3 




e-¥) 








~ 2 


Anma 



CA 



3n 2 m 



(2d) 
(3d) 



(28) 
(29) 



with the cutoff energy u>a — — A 2 /m and the two-particle binding energy eg. 
Then the stress correlation function including the cutoff term in Eq. ([2~T|) reads 



Xxyxy (0) — P 
Xxyxy(0) = S — 



3C 
8nm 
C 
Auraa 



C UJ A 
In — 

47rm sb 

AC y/mujA 
157r 2 m 



(2d) 
(3d). 



(30) 
(31) 



The zero-range interaction leads to universal high-frequency tails rj(uj) ~ C/(8raui) in 
2d [T71H5] and 7y(w) ~ C /(l&ir^Jmu)) in 3d [TBIUS] . These tails have to be subtracted 
to make the frequency integral (|22p convergent, and one obtains the shear viscosity 
sum rules UMIMIHI 



dui 
dui 



r)(uj) - 



C 
8raui 
C 



0(oj-e B ) 



= P- 



3C 
8ura 



= e — 



C 
8ura 



I5iry/mu! 



= P- 



C 



_ 2 
Airma 3 



C 



6nma 



(2d) 
(3d). 



(32) 
(33) 
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The universal high-frequency behavior is most clearly seen if one looks near the quan- 
tum critical point at zero density and zero temperature [241125] . The shear viscosity 
in this limit but with the same value of C as in the dense system has the form [19] 

By subtracting ?7o(w) one arrives a low-energy sum rule which captures only the 
finite-density effects 19 

2 f°° 

- / dw[r)(uj)-r)o(u)] = P (2d). (35) 

T JO 

In conclusion, we have argued that zero-range interactions realized in ultracold 
atomic systems do not modify the spin /-sum rule but lead to characteristic contact 
terms in the shear viscosity. 
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